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Abstract 

The reaction pp — » K + + (Ap) was measured at T p = 1.953 GeV and = 0° with 
a high missing mass resolution in order to study the Ap final state interaction. The 
large final state enhancement near the Ap threshold can be described using the stan- 
dard Jost-function approach. The singlet and triplet scattering lengths and effective 
ranges are deduced by fitting simultaneously the Ap invariant mass spectrum and 
the total cross section data of the free Ap scattering. 
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1 Introduction 

The study of the Ap interaction is part of the systematic investigation of the 
hyperon-nucleon (YN) interaction. The YN interaction is an ideal testing 
ground for the role of strangeness in low and intermediate energy physics. It is 
also of great relevance in studies of the SU(3)fl av0 r symmetry. In addition the 
YN interaction is needed for studies of hypernuclei. Experimental information 
on the Ap interaction has been obtained from Ap scattering experiments [Tf2"] , 
the binding energies of hypernuclei [3"]T4"] . and studies of the Ap final state in- 
teraction (FSI). The Ap FSI was observed in strangeness transfer reactions 
K~ + d — > 7r~ + (Ap) [5] and in reactions with associated strangeness produc- 
tion 7T+ + d -> K+ + (Ap) and p + p ->■ K+ + (Ap) PMM111T2TTM1IT5] . The 
first high resolution study of the reaction p + p — » K + + (Ap) was performed 
at Saclay [H] with proton kinetic energies of 2.3 and 2.7 GeV and missing 
mass resolutions between 2 and 5 MeV depending on the scattering angle. 
The FSI enhancement close to the Ap threshold was later analyzed [16] in 
terms of the inverse Jost function and the effective range approximation |17j . 
Recently a new method based on dispersion theory was developed [T8|ri9] in 
order to extract the scattering length and the effective range and to estimate 
the theoretical error. 

Theoretical studies of the Ap interaction with the Nijmegen [20] and Julich 
[21~] meson-exchange models predict the potentials, phase shift parameters 
and effective range parameters. Recently, the Ap interaction has been studied 
using the chiral effective field theory [221123] . Another topic of the Ap FSI is 
the prediction of a narrow S=-l dibaryon resonance by the cloudy bag model 
with an invariant mass of 2109 MeV [24] . 



2 Experiment 

The reaction p + p — > K + + (Ap) was measured at 0° using the proton beam 
from the the cooler synchrotron COSY, the magnetic spectrograph BIG KARL 
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Fig. 1. Left: Layout of the magnetic spectrograph BIG KARL. The charged particle 
tracks are measured in the focal plane using two stacks of multiwire drift chambers, 
two threshold Cherenkov detectors and two scintillator hodoscopes. Right: TOF 
spectrum with pion suppression by two Cherenkov detectors. 

[25] and a 1.0 cm thick liquid hydrogen target (see Fig. [1]). Since the most im- 
portant part of the FSI enhancement is located in the sharply rising part near 
the Ap threshold a high missing-mass resolution was required. The momen- 
tum of the incoming proton beam was about 2.735 GeV/c corresponding to a 
kinetic energy of 1.953 GeV. The absolute beam momentum was found from 
the fits of the kinematic parameters for two reactions: p + p — > K + + (Ap) 
and p + p — > d + tt + where kaons and deuterons were measured simultaneously 
at BIG KARL momentum 1070 MeV/c. The absolute precision of the beam 
momentum was 0.15 MeV/c. The scattered particles in the momentum range 
930 - 1110 MeV/c were detected in the focal plane using two stacks of mul- 
tiwire drift chambers, two threshold Cherenkov detectors and two scintillator 
hodoscopes. The ratio of beam momentum to scattered particle momentum 
was ideally suited for a measurement at 0°. In the first dipole magnet of BIG 
KARL the beam was magnetically separated from the scattered particles and 
guided through the side exit of the outer yoke into the beam dump. Thus, the 
huge background from dumping the beam within the spectrometer has been 
avoided. Particle identification was performed using the energy loss (Ai?) and 
time of flight (TOF) information from the scintillator hodoscopes. In addi- 
tion two threshold Cherenkov detectors [26] were used in order to achieve a 
very high pion suppression factor of 10 5 . The momentum of the kaon was 
measured and the missing mass of the Ap system was deduced. The effective 
missing mass resolution depending on the spread of the beam momentum, 
the momentum resolution of the magnetic spectrograph and the 1 MeV bin 
width was a = 0.84 MeV. In order to cover the missing mass range 2050 - 
2110 MeV three overlapping settings of the spectrograph (mean momenta: 
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1070, 1010 and 960 MeV/c) were used. The relative precision of the momenta 
of the three settings was 0.1 MeV/c. 



Acceptance corrections with respect to solid angle and momentum were taken 
from Monte Carlo calculations. The acceptance correction functions contained 
the magnetic spectrograph momentum acceptance around 0° with (p x ,P y ,Pz) 
cuts which for a given missing mass bin corresponded to the measured solid 
angle dQ. They contained also the detector efficiency corrections. The de- 
tector efficiency included efficiency of scintillator detectors which determined 
trigger and particle identification and magnetic spectrograph efficiency with 
field inhomogeneity at the edges of the acceptance. Acceptance corrections as 
determined by Monte Carlo simulations were checked by use of experimental 
distributions of simultaneously measured pions. They are shown in Fig. |2j The 
acceptance correction function of the 960 MeV/c setting looks different due 
to a slightly different (p x ,P y ,Pz) cut. 
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Fig. 2. Left: Acceptance correction functions. Right: Kaon survival probability 



The kaon decay along flight path was taken into account for each individual 
trajectory. Path lengths obtained from calculations of tracks in the magnetic 
field with the Turtle code (27] were checked by experimentally deduced values 
from time of flight measurements of pions and protons. The cross section error 
due to survival probability is less than 1 %. The final kaon survival probability 
averaged for a given missing mass bin is presented on the right side of Fig. |5J 

The relative normalizations of the three different spectrograph settings were 
deduced from luminosity monitors located in the target area which were inde- 
pendent of the spectrograph settings. The relative normalization errors due to 
the luminosity measurement were negligibly small. The relative normalization 
errors due to the acceptance corrections were estimated to be less than 2 %. 

The absolute cross section normalization was determined by measuring the 
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luminosity as described in [25] . At the beginning of each beam period we cali- 
brated two luminosity monitors counting the left and right scattered particles 
from the target as function of the number of beam particles. To this end, the 
beam current was highly reduced in order to count the beam particles individ- 
ually with a fast scintillator hodoscope in the beam dump. For the calibration, 
the dependence of the luminosity signal on the beam intensity was fitted by 
a linear function. The resulting beam intensity error amounted to 5 %. The 
density of the bubble-free liquid hydrogen target (p = 0.0775 g/cm 3 ) with 
1 /im thick mylar foil windows was kept constant by stabilizing the temper- 
ature to 15.0 ± 0.5 K using a high-precision temperature control [28]. The 
target thickness, i.e. the distance between entrance and exit window (nominal 
1 cm) was precisely measured with a calibrated optical telescope. The target 
thickness error was about 5 %. The overall systematic normalization error was 
estimated to be 10 %. The missing mass spectrum is shown in Fig. [3J 

q (MeV/c) 




M x (MeV) q (MeV/c) 

Fig. 3. Left: Missing mass spectrum of the reaction p + p — > K + + (Ap) measured at 
T p = 1.953 GeV and @k = 0°. The upper axis indicates the cm. momentum q of the 
Ap system. Solid line: Combined six-parameter fit. Dashed line: p + p — > K + + (Ap) 
phase space distribution. The region above 2090 MeV was measured with very small 
statistical errors in order to study a possible resonance anomaly near 2096.5 and/or 
2098.0 MeV [llj. Right: Total Ap — > Ap cross section |l|2j vs. cm. momentum q. 
Solid line: Combined six-parameter fit. Dotted line: Spin-averaged parameters. 



3 Formalism 

The observed missing mass spectrum can be described by factorizing the reac- 
tion amplitude in terms of a production amplitude and a final state enhance- 
ment factor. The method of parametrizing the FSI enhancement factor by 
the inverse Jost function [TF] is described in [IB] . Taking the spin statistical 
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weights into account the double differential cross section may be written as 
d 2 a 



dn K dM, 



Ap 



4 q 2 + aj 4 



q 2 + A 2 
q 2 + a 2 



(3.1) 



Here, |M S | 2 and \M t \ 2 are the singlet and triplet production matrix elements 
squared, q the internal c.m.-momentum of the Ap system, a s , (3 S , a t , (3 t the sin- 
glet and triplet potential parameters, and $3 the ratio of the three-body phase 
space distributiorJZI and the incident flux factor. The potential parameters a 
and (3 can be used to establish phase-equivalent Bargmann potentials [29]. 
They are related to the scattering lengths a, and effective ranges r of the low- 
energy S-wave scattering, a = (1 — \Jl — 2r/a) /r, j3 = (1 + yd — 2r/a)/r. 
Instead of the parameters a s , (3 S , a t and (3 t one can equally well use the sin- 
glet and triplet scattering length and effective range parameters a s , r s , a t and 
r t . The expression (13.11) is folded with the Gaussian missing-mass resolution 
function (a = 0.84 MeV) before comparing with the data in the fit program. 

The total cross section a of the free Ap elastic scattering can be expressed at 
low energies as a function of the c.m.-momentum q using the effective range 
approximation [IT] . 




4 Fit Results 



The nonlinear least-square fits are performed using the program Minuit from 
the CERN library. We determine the spin-averaged scattering length a and 
effective range f in a three-parameter fit by applying the constraints |M S | 2 = 
|M 4 | 2 =|M| 2 , a s —at—a, r s =r t =f. We fit only the missing mass spectrum with- 
out taking the Ap total cross section data into account. The fit yields an 
excellent description of the missing mass spectrum with xi-ed = 0.55, \M\ 2 = 
27.8j£g b/sr, a = -2.43ig;$ fm and f = 2.21±g;l| fm. But these spin-averaged 
parameters fail completely to reproduce the total Ap— >Ap cross section data, 
see dotted curve in Fig. [3] on the right (Xred = ^)- This result indicates that 
the singlet and triplet effective range parameters are different. 

Therefore, we study the missing-mass spectrum and the total Ap cross sections 
simultaneously in a combined fit. In addition we take a previous measurement 
of the Ap FSI in the reaction K~ + d — > 7r~+p + Aat rest into account 

1 Note the difference of factor 2 in the definition of $3 with respect to |30| 
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[5] which yielded an independent determination of the triplet parameters, 
a t = —2.0 ±0.5 fm and r t = 3.0 ± 1.0 fm. These two values and their errors are 
taken as additional experimental data, i.e. as 1-a constraints for the fitted a t - 
and r r values in the combined fit. In order to explore the parameter space we 
start with a five-parameter fit keeping the ratio \M t /M s \ 2 fixed. Starting with 
\M t /M s \ 2 = 1, we varied the ratio \M t /M s \ 2 between and 8. The resulting 
X 2 values depend rather strongly on the ratio \M t /M s \ 2 . The optimum fit is 
achieved with \M t /M s \ 2 = 0. A six-parameter fit confirms this result. The final 
results are listed in Table [TJ 

Table 1 



Six-parameter fit results. 



|M S | 2 (b/sr) 


a s (fm) r s (fm) 


\M t \ 2 (b/sr) 


a t (fm) 


n (fm) 


Xred 


mil* 


9 ,o+0.16 99-1+0.16 
z - 4O -0.25 z - zi -0.36 


n n+ 19 


l.OO_ 2 2 


o 7+O.6 
-J-' -0.6 


0.53 



We note that a possible theoretical uncertainty of the Jost-function approach 
may be in the order of 0.4 fm for the scattering length and even more for the 
effective range as suggested by the analysis of pseudodata [15] . This aspect 
will be discussed in a further paper. 

Concerning the experimental uncertainties, the overall normalization error 
yields a systematic error of 10 % for the production matrix element (|M| 2 
and |M S | 2 ). The effective range parameters are not affected. They depend only 
on the relative shape of the missing mass spectrum. The scattering length a 
is mainly determined by the strongly rising part of the spectrum near the 
Ap threshold. The effective range parameter r depends mainly on the spectral 
distribution towards higher invariant masses. The relative normalization errors 
between the three parts of the spectrum yield a systematic error of 0.02 fm for 
the scattering length (a and a s ) and 0.13 fm for the effective range parameter 
(f and r s ). These errors are smaller than the error estimates of the fit, but they 
must be taken into account in the evaluation of the total error. The precision 
of the beam momentum and the kaon momenta is so high that the deduced 
FSI parameters are not affected. 

The results shown in Table 1 depend on the accuracy of the included Ap data, 
especially on the overall normalization of the cross sections, see the right side of 
Fig. [31 The error bars indicate the statistical errors which are the main source 
of errors in those hydrogen bubble chamber experiments. The systematic errors 
are small compared to the statistical errors. The data are based on 378 and 224 
elastic Ap scattering events, respectively. They were taken by two independent 
groups [Tf2] using the 81-cm hydrogen bubble chamber at CERN. The two data 
sets are consistent within the errors. 
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5 Discussion and conclusion 



The reaction p + p — > K + + (Ap) was measured at T p = 1.953 GeV and 
O = 0° with a high missing mass resolution in order to study the Ap FSI. 
A three-parameter fit with spin-averaged effective range parameters a and f 
yields a good description of the missing mass spectrum but fails completely to 
describe the momentum dependence of the free Ap scattering. The combined 
study of the missing mass spectrum and the free Ap scattering using five- 
and six-parameter fits reveals that the production of the Ap system in the 
triplet state is rather small. The best fit yields \M t \ 2 = and the la limit 
is reached for \M t /M s \ 2 = 0.168 corresponding to |M S | 2 = 73.6 b/sr and 
\M t \ 2 = 12 A b/sr. In the fit we take also the independent determination of 
the triplet parameters a t and r t by Tai Ho Tan [5] into account. 

It is interesting to note that the spin-averaged parameters a and f deduced 
from a three-parameter fit of the missing mass spectrum are identical with the 
singlet parameters a s and r s deduced from a six-parameter fit of both the miss- 
ing mass spectrum and the free Ap scattering. Also, the fit indicates that the 
reaction p+p — > K + + (Ap) is dominated by the singlet contribution, however, 
this result should be considered with respect to the theoretical uncertainties 
of the Jost-function approach [19]. A direct determination of the singlet and 
triplet contributions requires polarization measurements as proposed in |18j . 

The result \a s \ > \a t \ indicates that the Ap interaction is more attractive in 
the singlet state than in the triplet state. This finding is in accordance with 
the analysis of the binding energies of light hypernuclei jl] . The present results 
agree within errors with recent predictions of the new Julich meson-exchange 
model J04c [21] yielding a s = —2.66 fm, r s = 2.67 fm, a t = —1.57 fm and 
r t = 3.08 fm. Recent calculations of the NLO effective field theory yield similar 
results for a s , r s , a t and r t [23J. The present results are also in agreement with 
a similar analysis of the reaction p + p — > K + + (Ap) at T p = 2300 MeV 
and G = 10.3° |16j . Thus, the method of factorizing the reaction amplitude 
in terms of a production- and FSI-amplitude yields comparable results for 
different bombarding energies and scattering angles. In this sense, the present 
study provides also a systematic check of the method. 

We analyze our data using the standard Jost-function approach [T7]. The 
recently proposed dispersion-integral method [T8iri9"] differs from the Jost- 
function approach by restricting the dispersion integral over the scattering 
phase shift to a finite upper limit. In a further investigation we will analyze 
our data using different approaches. 

A narrow S=-l resonance predicted by the cloudy bag model [23] is not visible 
in our missing mass spectrum. The small structures at 2096. 5±1. 5 and 2098. 0± 
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1.5 MeV observed in a previous experiment [TT] are not confirmed. Further 
data analysis is underway in order to quantify the resonance limits. But this 
is beyond the scope of the present paper. 
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